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It is shown that the method of infinitesimal generators ("Racah's method") can 
be broadly and systematically formulated as a method applicable to the calculation 
of reduced coupling coefficients for a generic subalgebra chain G D H, provided the 
' reduced matrix elements of the generators of G and the recoupling coefficients of 

H are known. The calculation of SO (5) D SO (4) reduced coupling coefficients is 
considered as an example, and a procedure for transformation of reduced coupling 
Q^, coefficients between canonical and physical subalegebra chains is presented. The 

fH I problem of calculating coupling coefficients for generic irreps of S0(5), reduced with 

I respect to any of its subalgebra chains, is completely resolved by this approach. 

lMj; PACS numbers: 02. 20.Qs 

^ ■ I. INTRODUCTION 

i> : 

^ i Continuous symmetries and their associated Lie algebras facilitate the description of 

^ ' many-body systems both directly and indirectly. When a symmetry occurs as a dynamical 
Q ■ symmetry of the system, the corresponding algebra immediately gives the spectroscopic 
O . properties of the system. However, even when a symmetry is strongly broken, the algebraic 
I structure nonetheless provides a calculational tool, classifying the basis states used in a full 
. ^ ■ computational treatment of the many-body problem and greatly simplifying the underlying 

^ ' calculational machinery. Lie algebras have a long history of application, in both these 
^ . capacities, to nuclear spectroscopy and related problems. The fundamental quantities 
underlying calculations within a Lie algebraic framework are the coupling coefficients of the 
algebra, also known as generalized Clebsch-Gordan coefficients or Wigner coefficients. These 
are needed in order to couple states (or operators) of good symmetry to yield new states (or 
operators) of good symmetry, and they are required for the calculation of matrix elements 
through the generalized Wigner-Eckart theorem of the algebra. 

The Lie algebra SO (5), isomorphic to Sp(4), has several distinct applications in nu- 
clear theory, involving different physical realizations of the operators, and in which different 
subalgebra chains are relevant to the symmetry properties. The natural construction of 
SO (5) in terms of generators of rotation in five- dimensional space gives rise to a canonical 
S0(4) ~ S0(3) ® S0(3) subalgebra."^ However, application as the proton-neutron pairing 
quasispin algebra^"^ requires reduction with respect to the U(l) ® S0(3) algebra of isospin 
and occupation number operators. For the dynamics of spin-2 bosons (as in the interacting 
boson model^°'^^) and for the Bohr collective model, ^^"^^ the appropriate reduction is instead 
with respect to a physical angular momentum SO (3) subalgebra. 
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In this article, it is shown that the method of infinitesimal generators ("Racah's method") 
can be broadly and systematically formulated as a method applicable to the calculation of 
reduced coupling coefficients for a generic subalgebra chain G D H, provided the reduced 
matrix elements of the generators of G and the recoupling coefficients of H are known 
(Sec. II). More specifically, the problem of calculating coupling coefficients for generic irreps 
of SO (5), reduced with respect to any of the subalgebra chains, is completely resolved by 
this approach. The calculation of reduced coupling coefficients for the SO (5) D SO (4) 
canonical chain is considered in detail (Sec. III). Couphng coefficients reduced with respect 
to the noncanonical subalgebra chains of S0(5) may be obtained by a similar application of 
Racah's method, or they can be deduced from the canonical chain coupling coefficients by 
unitary transformation. The general formulation in the presence of outer multiplicities for 
H, numerical examples for SO (5), and a detailed account of the transformation procedure 
between subalgebra chains are given in the appendices. 



II. METHOD 

A. Background and definitions 

Consider a Lie algebra G and subalgebra H. States which reduce this subalgebra chain 
may be identified by the irrep labels F of G, the irrep labels A of i7, and a label A (typically 
the Cartan weights) to distinguish basis states within A, as | FAA). The coupling coefficients, 
or generalized Clebsch-Gordan coefficients, for G relate the uncoupled product states of two 
irreps of G to the coupled states, as 



^ A^ A^ 
A1A2 V'^i ^2 

A1A2 




(1) 



In general, additional labels will be required to resolve multiplicities. There may be "outer" 
multiplicities in the Clebsch-Gordan series for the outer product of G (that is, Fi ® F2 may 
contain the irrep F more than once), and there may be "branching" multiplicites under the 
restriction oi G to H (that is, the given irrep T oi G may contain an irrep A of H more than 
once). The coupling relation (1) generalizes, with multiplicities, to 
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aiAi 02^2 
Ai A2 



(2) 



where p is the outer multiplicity index for G ® G — )■ G, and the a indices resolve the 
branching multiplicities for G ^ H. Furthermore, H may be subject to outer multiplicites 
[H ®H ^ H). In the following discussion, we shall for simplicity take the subalgebra H to 
be multiplicity free. Such is the case for the commonly encountered situation in which the 
physically relevant subalgebra H is SO (3), as well as for the subalgebra SO (4) considered 
in Sec. III. However, the necessary generalizations in the presence of outer multiplicities on 
H are given in Appendix A, as would be needed for consideration of, e.g., chains involving 
SU(3) as a subalgebra. 
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Racah's factorization lemma^^ allows the coupling coefficient appearing in (2) to be de- 
composed as the product 



Ai A2 



Ai A2 
Ai A2 



A\ / Ti 12 
A / 1 aiAi a2A2 



aA 



(3) 



of a coupling coefficient of H, embodying all the dependence upon weights A, with a reduced 
coupling coefficient (or isoscalar factor) for G D H. The reduced coupling coefficient is 
nonzero only if F is contained in the outer product of Fi and F2 (z.e., Fi (8) F2 — >■ F), each 
irrep of H is contained in the corresponding irrep of G {i.e., Fi Ai, T2 A2, and F — )■ A), 
and A is contained in the outer product of Ai and A2 {i.e., Ai (8) A2 ^ A). The reduced 
coupling coefficients satisfy the orthonormality conditions^® 
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and 



V( Ti r2 
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(5) 



for any irrep A such that F — > A. 

If T^^ is an irreducible tensor operator with respect to H, the Wigner-Eckart theorem 
for H permits the expression of a general matrix element of T^^ as^® 
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A At 
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(6) 



in terms of a coupling coefficient for H and a reduced matrix element with respect to H. A 
Wigner-Eckart theorem of this form [or its generalization (A4)] may be obtained whenever 

H is a compact, semi- simple Lie algebra. 

Several methods may be considered, in general, for constructing the reduced coupling 
coefficients of Lie algebras: 



(1) Recurrence relations among coupling coefficients may be obtained by considering the 
action of an infinitesimal generator Gi = G[^'' + G^^ on uncoupled and coupled states. 
This approach, used in the present construction, is broadly termed "Racah's method" 
(see Ref. 18) and generalizes the classic recurrence method for evaluating SU(2) ~ 
S0(3) Clebsch-Gordan coefficients.^^ 

(2) Recurrence relations and seed values may be obtained by considering the action of 
a "shift tensor", lying outside the algebra, which connects different irreps of the 
algebra.^° 

(3) Consistency relations among coupling and recoupling coefficients serve as the basis for 
a "building up" process, ^'^^ in which unknown coupling coefficients can be deduced 
from a few known coefficients. 
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(4) Explicit realizations of an algebra can be obtained in terms of bosonic or fermionic cre- 
ation and annihilation operators. Relations among coupling coefficients follow from 
considering the matrix elements of tensor operators acting on bosonic or fermionic 
states {e.g., Ref. 22). This approach is generally restricted to symmetric irreps, anti- 
symmetric irreps, or irreps which can be obtained as simple combinations thereof. 

Indeed, all of these approaches have been applied or suggested, in various forms, for the 
calculation of specific classes of S0(5) coupling coefficients. ^'^'^^'^^"^^ For the symmetric 
irreps of SO (5), one may also work with an explicit realization in terms of five-dimensional 
spherical harmonics as functions on the four-sphere. Their triple overlap integrals are then 
proportional to SO (5) coupling coefficients. ^*^'^° 



B. Racah relations among reduced coupling coefficients 



Let us now consider how the first approach, i.e., Racah's method^'^^ based on the action 

of infinitesimal generators, can be generally and systematically formulated as a method 
applicable to the calculation of reduced coupling coefficients involving generic irreps of an 
arbitrary subalgebra chain. Consider the action of an infinitesimal generator Gi of G on 
the coupled product state of (2). The generator on the product space is of the form Gi — 
G\^^ + Gf\ where Gf^ acts only on the space carrying the irrep Fi and acts only on 
the space carrying the irrep The equivalence of the action of Gi on the two sides of (2) 
imposes conditions on the coupling coefficients connecting the different basis states used on 
the two sides. For effective application of Racah's method, it is most convenient to recast 
these relations among coupling coefficients so that they involve only (1) reduced coupling 
coefficients of G with respect to H, (2) reduced matrix elements of the generators of G, and 
(3) recoupling coefficients of H, as obtained in this section. 

Racah's approach requires that the action of the generators on the basis states of an irrep 
be known explicitly. In general, if coupling coefficients are to be determined for states which 
reduce G D H, it is necessary to consider the action of the generators which are in G but 
not in H, since only these generators can connect different irreps A of H. Note that Racah's 
method is essentially an extension of the classic scheme^^ for calculating the ordinary SO (3) 
Clebsch-Gordan coefficients, via recurrence relations obtained by considering the known 
matrix elements of J± — -|- j]f ^ , between the uncoupled product states | ) and 
coupled states | ^) (see also Sec. II C). 

If T^J is a generator of G, expressed as an irreducible tensor operator with respect to H, 



we begin by considcrin 
and coupled product states. 



the matrix element of T^'^ 

At 



j^At (1) 
At 



At (2) 
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Ai A2 



At (2) 




(7) 



The coupling relation (2) may be used to express each ket on the right hand side of (7) 
entirely in terms of uncoupled states, and its inverse, obtained by orthonormality of coupling 
coefficients, may be used to express the bra on the left hand side entirely in terms of coupled 
states. Since T^'^, as a generator of G, does not connect different irreps of G, and since the 
matrix elements of T^'^ between states within an irrep of G depends only upon the irrep 
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labels, the result simplifies to 
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To introduce reduced coupling coefficients and reduced matrix elements, we apply Racah's 
factorization lemma (3) and the Wigner-Eckart theorem (6), yielding 
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It now remains to ehminate the coupling coefficients of ff, and thus all reference to weights. 

The orthogonality relations allow these coefficients to be moved to the right hand side, 
resulting in sums of quadruple products of coupling coefficients. These sums arc recognized 
as recoupling coefficients of if, specifically, the "unitary 6-A symbols", or transformation 



brackets between basis states in the couphng schemes [(AiA2)^^^A3 
which are given by 



and [Ai(A2A3)^-]^, 
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The 6-A symbol is nonvanishing only if the Clebsch-Gordan series relations Ai A2 ^ A12, 
A12 ® A3 A, A2 ® A3 A23, and Ai (g) A23 A are satisfied. Let $(AiA2; A) denote the 
phase factor incurred by interchange of the first and second irreps in a coupling coefficient 
of H, i.e., {^^ Ai I a) = *(A2Ai; A)(^J A2 | a)- Then the condition (9) becomes, with labels 
renamed for simplicity. 
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expressed entirely in terms of the reduced coupling coefficients to be calculated, reduced 
matrix elements, and recoupling coefficients of the lower algebra H. 

For the important special case in which H is the angular momentum algebra SO (3) ~ 
SU(2), the relation (11) becomes 



E 



r 

aJ 



r \ T, T, 



aJ 



E( \J2+Jt-\-J' J Jl] J2 J'l J'\ Ti 
\Jt J Ji\\aiJi 



a[J[ 
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where J = (2J + 1)^/^. Note that the customary form^^ of the Wigner-Eckhart theorem 
for S0(3) is defined in terms of a 3-J symbol rather than a Clcbsch-Gordan coefficient. 
This results in a reduced matrix element which differs in normalization and phase from the 
definition implied by the generic statement of the Wigner-Eckhart theorem in (6). The 
reduced matrix elements under these two conventions are related by (J3 || T^-^^) || J-^^q^^,) — 

{-f''MM\T^''^\Ji)- 



C. Solution of the homogeneous system 

The condition (11) yields a different relation among specific reduced coupling coefficients 
for each choice of values for the four irrep labels aiAi, a2^.2-i A, and a' A' (the multiphcity 
index a is summed over). If there are A'^ coupling coefficients for the coupling Fi (g) r2 — > pF, 
then the relations obtained from (11) constitute a linear, homogeneous system of equations 
in A^ unknowns for these coupling coefficients. 

Note that the most familiar and traditional approach to extracting coupling coefficients, 
after obtaining some set of relations among them, is to proceed by recurrence {e.g., in the 
familiar case of S0(3)^^ and in the "building-up process", as well as in prior applications 
of Racah's method to higher algebras^^'^^). That is, a seed value is given for one coupling 
coefficient, and further coefficients are deduced inductively (one by one) from those already 
obtained. 

A recurrence approach is indeed natural in the case of S0(3). The relations obtained by 
considering the actions of J± are 

K^{JM) ( I ^J^, ) = X±( JiMi T 1) ( M%i A Uf ) + K^{J2M2 T 1) ( ii, m.^i U ) , (13) 

where M = M1+M2TI, and KJ JM) = {JM ± 1 1 J± | JM) = [( J^M)( J±Af +1)]!/^ is the 
generator matrix element. These relations connect at most three coupling coefficients, and a 
natural order for traversing the coefficients can easily be chosen, such that only one unknown 
arises at each step, as illustrated in Fig. 1(a). [This is accomphshed by involving certain 
known-zero or "forbidden" Clebsch-Gordan coefficients, represented by the triangle vertices 
without dots in Fig. 1(a), in the relations.] Since classic treatments {e.g., Ref. 31) apply 
orthonormality relations interspersed with the recurrence relations at intermediate stages of 
the calculation, we stress that recourse to orthonormality conditions is not actually necessary. 
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FIG. 1: The classic problem of constructing the S0(3) Clebsch-Gordan coefficients ( | M1+M2 ) 
by use of the relations (13). Dots indicate allowed non-zero coefficients. Coefficients at the vertices 
of a dashed triangle are connected by (13). (a) The conventional recurrence approach, in which 
coefficients are calculated inductively from a seed coefficient, making use of relations which in some 
cases invoke known-zero Clebsch-Gordan coefficients, (b) A full set of relations among allowed 
Clebsch-Gordan coefficients, yielding a linear, homogeneous system of equations in the Clebsch- 
Gordan coefficients. 



It may be seen from the figure that all coefficients are accessible by the relations (13). In 
anticipation of the treatment of higher algebras, we also observe that all allowed coefficients 
may be connected by the relations directly, without involving any forbidden coefficients, as 
in Fig. 1(b). This yields a system of equations which fully determines the coefficients, to 
within an overall phase and normalization, although in this case the system is not amenable 
to solution by recursive calculation of successive coefficients from a single seed coefficient. 

For higher algebras, many irreps of H may be connected by the generator T*^^^\ and 
therefore each relation obtained from (11) may involve many unknown coupling coefficients. 
A simple recurrence pattern, as in Fig. 1(a), may be impractical to devise. A more generally 
applicable and straightforward approach is to directly solve the linear, homogeneous system 
of equations for the unknown coupling coefficients, by standard linear algebraic methods, 
e.g., Euler row reduction. 

Let us therefore summarize the system of equations which must be constructed and 
solved. The N unknown couphng coefficients for Fi (g) r2 pT may be labeled with a 
single counting index as = {f^\^ ^^^2 I oa) ^^^^ * ~ 1, . . . ,N. Each value of the index i 
therefore designates a specific combination (aiAia2A2oA). The numerical coefficients of the 
unknown quantities Cj in the relations (11) do not depend upon the outer multiplicity index 
p, to be discussed further below. For each generator T^^"^^ in G but not in H, and for each 
quadruplet of irrep labels (aiAia2A2Aa'A'), the relation (11) yields an equation (which we 
label by a counting index k) of the form Xli^i '^kiCi = 0, i.e., linear and homogeneous in the 
Ci. The equation is nonnull {i.e., some of the coefficients a^i are nonvanishing) only if the 
Clebsch-Gordan series conditions Ai ® A2 — )■ A and A' ® A^ A are met.'^^ 

The resulting equations must be aggregated to yield the full system, which may be ex- 
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pressed in matrix form as 



Coupling coefficient 

(o]^A]^a2A2aA) 




roi 







(14) 



=A 



Normally, it suffices to consider the conditions obtained with (A1A2AA') such that Fi 
Ai, r2 A2, r — )■ A, and F — > A', that is, relations involving only "allowed" coupling 
coefficients, as in Fig. 1(b). However, in certain exceptional cases, additional conditions 
involving known-zero coupling coefficients may be necessary, analogous to Fig. 1(a). These 
may be obtained, e.g., by considering some A' with F ^ A'. 

The problem of solving this linear homogeneous system of equations (14) is equivalent 
to finding the null vector (or vectors) of the matrix A appearing on the left hand side 
of (14). In general, there may be many more rows (equations) than columns (unknown 
coupling coefficients). However, these rows are not linearly independent. In the case where 
Fi ® F2 — )■ F is free of outer multiplicity, the matrix A can be expected to be of rank 
N — 1. The null vector is then uniquely determined, to within normalization and phase, and 
its entries arc the coupling coefficients [C1C2 ■ ■■ C^]. The proper normalization, yielding 
coefficients satisfying the condition (4), is obtained by evaluating 



That is, the summation runs over the subset of entries Q sharing the same value for aA. An 
identical result for J\f must be obtained, regardless of the choice of aA, provided F — > oA. 
(In fact, the requirement of equality may be used as an internal consistency check on the 
calculation.) The normalized coupling coefficients are then obtained by dividing the null 
vector by J\f. 

An overall sign remains to be chosen for the entire set of coupling coefficients for Fi(8>F2 
F. For instance, Refs. 23,24 suggest a"generalized Condon-Shortley phase convention", such 
that {I'^fjJD > 0' that is, a positive value is adopted for the coupling coefficient 
involving the highest weight irreps of H contained in Fi and F and the highest weight irrep 
A2 consistent with these. 

More generally, when the couphng Fi (g) F2 — >■ F has outer multiphcity D (p = 1, . . . , D), 
the matrix A may be expected to be of rank N — D. That is, the system of equations 
given by (14) yields D linearly independent null vectors (or its null space has dimension D). 
The null vectors obtained by Euler row reduction must be orthonormalized,^^ e.g., by the 
Gram-Schmidt procedure, to yield a set of coupling coefficients satisfying the orthonormality 
relation (4). Note that the appropriate inner product for this orthonormalization is not the 
standard vector dot product on R^. Rather, if we label the entries of each null vector Cp 
as [Cpi Cp2 • • • C'pjv] , then the inner product to be used for orthonormalization is 




(15) 



(same aA) 




(16) 



(same oA) 
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The same value of M.pip is obtained regardless of the choice of oA used in evaluating the 
sum. (Again, requirement of this equality provides an internal consistency check on the 
calculation.) The orthonormal coupling coefficients are then simply the entries of the or- 
thonormalized null vectors. 

When an outer multiplicity is present, it should be noted that the coupling coefficients are 
defined only to within a unitary transformation, arising from the arbitrariness inherent in 
defining the resolution of the outer multiplicity, i.e., in choosing the basis states \ pT ■ ■ ■) (p = 
!,...,£>) spanning the D-dimensional space of irreps of type F. In the present calculational 
procedure, the freedom in resolution of the multiplicity is manifested in the freedom to 
choose different sets of orthogonal basis vectors for the null space of A. 

III. COUPLING COEFFICIENTS FOR SO(5) IN THE CANONICAL BASIS 

A. Overview 

For a concrete example of the apphcation of Racah's method in terms of reduced couphng 

coefficients, as developed in Sec. II, we consider the calculation of coupling coefficients for 
SO (5), reduced with respect to the canonical subalgebra SO (4). That is, we have SO (5) D 
SO (4) as the algebras G D H. Both essential criteria for application of the method are 
met: (1) the couphng and recoupling coefficients (Wigner calculus) for S0(4) are known, 
and (2) the reduced matrix elements of the S0(5) generators, considered as tensor operators 
with respect to S0(4), are also known. 

The algebra SO (5) contains several subalgebra chains, involving distinct SO (3) ~ SU(2) 
subalgebras, 

S0(5) D S0(4) D S0j(3) D S0j(2) 

[hh] [pq] J Mj 

D S0(4) ~ S0x(3) <8) S0y(3) D S0x(2) ® S0y(2) 

X Y Mx My 

D Ujv(l) (g) S0t(3) D S0t(2) 

K Ms T Mt 

D SO43) D SO42) 
a L Mi^ 

where the irrep label has been noted beneath each subalgebra. Branching multiplicity labels 
are indicated by k and a in the last two chains. Chain (I) is the standard canonical chain, 
while in (I') the canonical SO (4) subalgebra is reexprcssed using the isomorphism SO (4) ~ 
SO (3) ® SO (3). [As far as definition of reduced coupling coefficients is concerned, the 
two chains (I) and (F) are equivalent, but branching rules, coupling coefficients, etc., are 
simpler when expressed with respect to the latter chain (F).] The prerequisite definitions 
and algebraic results are summarized in Sec. IIIB, and the calculation of reduced coupling 
coefficients for the canonical chain is discussed in Sec. IIIC. 

Physical applications require the couphng coefficients of SO (5) reduced with respect to the 
noncanonical subalgebras of chains (II) and (III). The isospin algebra S0r(3) of chain (II) 
is the relevant subalgebra for the description of proton-neutron pairing. In this con- 
text, the S0(5) generators arise as quasispin operators for pairing of protons and neutrons 
occupying the same j-shell. On the other hand, the "physical" or "geometric" angular 
momentum subalgebra 80^(3) of chain (III) is the relevant subalgebra for apphcation to 



(I) 
(I') 
(11) 
(HI), 



(17) 
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systems of spin-2 bosons^°'^^ or the nuclear collective model. ''^^"''^^'^^ Explicit constructions 
of these subalgebras and further algebraic properties for the noncanonical chains are de- 
tailed in Appendix C, where the transformation between canonical and noncanonical bases 
is considered. 



B. Definitions and algebraic properties 

Let us begin with a concise but comprehensive summary of the construction of SO (5) 
and the algebraic properties needed for the application of Racah's method. Such a review 
is particularly necessary since notations and conventions for nearly all aspects of the treat- 
ment of S0(5) vary widely {e.g., Refs. 18,23,37,41), and phases and normalizations play an 
essential role in the calculation of coupling coefficients. 

The basic construction proceeds from the generators of rotation, 

L^g = —i(xrds — Xgdr). (18) 

These operators are Hermitian (Lj^ = L^s), are antisymmetric in the indices, and have 
commutators 

["^PQ'' -^rs] — ^(^^qr-^ps ~t~ ^ps-^qr ~t~ ^sq-^rp ~t~ ^rp-^sq) • ("^^) 

First, for S0(4), let Jr = \erstLst and Nr = Lr4 (1 < t, s,t < 3), i.e., 

Jl = -^23 J2 = ^31 = L12 , , 

Ni^Lu N2^L24 N^^L^i. ^ ' 

Then the Jr span the usual three-dimensional angular momentum algebra, which we de- 
note by S0j(3). The Jr and Nr together span S0(4), with commutators [Jr,Js] = iSrstJt, 
[Nr,Ns] = iErstJt, a^d [Jr, A^s] = iErstNt- The standard Cartan weight operators for S0(4) 
are J3 and A3. An SO (4) irrep is labeled by the highest weight defined by these operators, 
which is of the form [pg], with p > \q\, both integer or both odd half integer. 
The isomorphism SO (4) ~ S0x(3) S0y(3) is realized by taking 

Xk=l{Jk + Nk) Yk=\{Jk-Nk), (21) 

so [Xr, Xg] = iSrstXt, [Yr, Yg] = iSrstYt, and [Xr, Yg] = 0. The ladder operators for each S0(3) 
algebra are thus X± = Xi ± 1X2 and y± = Yi ± iY2. The natural Cartan weight operators 
in this scheme are then the S0(3) angular momentum projections Xq = X^ and Yq = Y3, 
defining weight labels Mx and My. An SO (4) irrep is then labeled by the highest weight 
(XY), i.e., the angular momenta associated with the S0x(3) and S0y(3) subalgebras. The 
S0x(3) ® S0y(3) irrep labels are related to the standard SO (4) labels hj X = ^{p + q) and 
Y = |(p — q) or, conversely, [p.q] = [X + Y,X — Y]. Note that the canonical S0j(3) is 
obtained as the sum angular momentum algebra of S0x(3) and S0y(3), since Jk = Xk + Yk, 
and thus the basis states reducing chains (I) and (I') are related to each other by ordinary 
angular momentum coupling. 

The Clebsch-Gordan series and coupling and recoupling coefficients for SO (4) follow im- 
mediately from the SO (3) ® SO (3) structure, most transparently with the (XY) labeling 
scheme for the irreps. The weights contained within (XY) are Mx — —X, . . . ,X — 1,X and 
My = — y, . . . ,Y — 1,Y. The Clebsch-Gordan series is given by apphcation of the triangle 
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inequality separately to each of the S0(3) algebras, that is, for (XiYi) ® (X2I2) (^^)) 
X= IX1-X2I, |Xi-X2| + l,...,Xi + X2 and y= 1^1-^2!, 1^1-^21 + 1,..., ^1 + ^2. Hence, 
no inner or outer multiplicities are obtained for SO (4). Coupling coefficients factorize into 
products of ordinary S0(3) Clebsch-Gordan coefficients, as 



(XiFi) (X2y2) 

MxiMyi Mx2My2 



(XY) 
MxMy 



Xi X2 
Mxi Mx2 



X 
Mx 



Yi ¥2 

Myi My2 



Y 

My 



(22) 



By inspection of (10), it is immediately apparent that the recoupling coefficients factorize 
as well, as 



iX^Y,) {X2Y2) {X,2Yu) 

[X-iY-i) (XY) (^23^23) 







X2 


Xl2 




Vi 


Y2 


Yi2 




Xs 


X 


X23 






Y 


Y23 



= (-) 



X1+X2+X3+XI 



fx. 


X2 


Xl2\ 




Y2 


YuX 


\X3 


X 


X23j 


IS 


Y 


Y23j 



(23) 



An equivalent result is given with standard S0(4) labels in Ref. 42. However, note that 
the result is considerably more cumbersome to derive if one uses the standard canonical 
chain (1).'^^ 

The algebra SO (5) includes the additional four generators (r = 1, . . . ,4). A tensor 
operator with respect to SO (4) is simply a simultaneous spherical tensor with respect to both 
the S0x(3) and S0y(3) algebras, i.e., a spherical "bitensor". For the S0(5) generators, we 
have bitensor expressions^^ 



-'^iio = X±i - T:^[{L23 



Lu) ± iiUi + L24)] 



>^o2? = ^±1 = Ti73[(^23 - L,,) ± t{L,, - L24)] 



(--) 



^2^2 
(--) 

rp(22' — rp 
2^2 



|(-^35 — i-^45) 



x^(lO) _ ^ 

- Y - 
^00 = -^0 - 

(--) 

+ 2 2 
(--) 

~2~2 



- \{Li2 + -^34) 
= \{Li2 — L^a) 

^l{Ls, + iL,,) (24) 

— |(-^15 ~ ■^-^25) • 



The phases are chosen so that these operators obey A^^^U^ 
generalization of the usual condition for a sclf-adjoint spherical tensor. All commutators 
involving X^ or have the values implied by the spherical bitensor notation of (24), e.g., 

[X±i,A 



'_\Mx+My a(^^) 



a 



(AV)i 



[|(A =F ;u)(A ± ;U + l)]^^^^(^±i)^/- The commutators between components of 
7^(1/21/2) given explicitly in Table I. 

The root vector diagram of S0(5) is shown for reference in Fig. 2(a), with the gener- 
ators (24) placed according to their S0x(3) ® S0y(3) weights [MxMy). The canonical 
subalgebra is highlighted in Fig. 2(b), and the construction of the physical subalgebras, of 
chains (I) and (HI), is indicated in Fig. 2(c,d) (see Appendix C). 

The standard Cartan highest weight labels for an SO (5) irrep are defined with respect to 
weight operators J3 and A^3 and have the form [/1/2], with li > I2, both integer or both odd 
half integer. It is more convenient in the present context to label SO (5) irreps by the highest 
weight defined by the S0x(3) (X>S0y(3) weight operators Xq and Yq (following Hecht^^). The 
resulting label has the form (RS), with R> S, each independently either integer or odd half 
integer. This label may also be considered as representing the angular momenta (X^Y^) 
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TABLE I: Commutation relations between the components of T^l'^^l'^ ^ for the S0(5) generator 
normaUzation and phase conventions defined in (24). 







7+- 




r__ 







~ 73^+1 


-73^+1 


-|(Xo + Fo) 









+3(^0-1^0) 








-i(Xo - Fo) 









+i(Xo + yo) 











FIG. 2: Root vector diagram for S0(5) and its subalgebras. (a) The generators of S0(5), labeled 
by their Cartan weights Mx and My. (b) The canonical subalgebra S0(4) ~ S0x(3) (g) S0y(3), 
which begins chains (I) and (I'), (c) The UAr(l) SOrCS) subalgebra of chain (II). (d) The 
physical angular momentum 80^(3) subalgebra of chain (III). These generators are obtained as 
linear combinations of the canonical generators, as indicated by the dashed lines, such that all 
generators have good S0l(3) D 80^(2) tensorial character. 
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TABLE II: Labeling schemes for irrcps of S0(5) in use in the physics hteraturc. with relations for 



intcrconvcrsion. 



Labels 



Range 



Relations 



Description 



Refs. 



[^1/2] 

(0102) 
(vf) 

m 
(a'ia'2) 

{RS) 



h— hih + 1, ■ 
h= 0, |, . . . 
ai= 0, 1, . . . 
a2= 0, 1, . . . 
v=0,l,... 
/=0, i,... 

^1= ^2' ^2 "t" -'-5 ' 

1^=0,1,... 
a'i= 0, 1, . . . 
4= 0, 1, . . . 
R= (S, + ^ , . 
S= 0,i,... 



ai= h — I2 
a2= 2I2 

v= h — I2 

f= h 

l[=h+l2 

^2= h ~ h 
a'i= 2I2 
a'2= h — I2 

R=^{h+l2) 
S=\{h-l2) 



S0(5) Cartan highest weight 6,37 

SO (5) Dynkin 46 

SO (5) Dynkin (modified) 47 

Sp(4) Cartan highest weight 5,45 

Sp(4) Dynkin 8,41 

S0(3) (8) S0(3) highest weight 23,48 



of the highest weight S0x(3) ® S0y(3) irrcp contained in the SO (5) irrep. The relation to 
the standard labels is i? = + ^2) and S = ^{li — I2) ■ A plethora of labeling schemes 
for SO (5) irreps are in use in the physics literature, interrelated as summarized in Table II 
(even more schemes arise if we consider the translation to physical labels, such as reduced 
isospin^^). 

The branching rule for S0(5) to S0(4), i.e., {RS) {XY), is given by23'37 



X ^R- 



n 



n 



(25) 



with < n < 2{R — S) and < m < 25", m and n integers. Graphically, the SO (4) irreps 
form a lattice bounded by a tilted rectangle, as illustrated by the dotted line in Fig. 3(c). 
The rectangle's "right" corner is at the highest weight (RS) = (Xj„y^), its "bottom" corner 
lies on the Mx axis, and the remaining two corners are specified by symmetry about the line 
Mx = My. The branching rule is shown for example irreps of SO (5) in Fig. 3: symmetric 
[Fig. 3(a)], antisymmetric [Fig. 3(b)], and generic [Fig. 3(c)]. 

The Clebsch-Gordan series for S0(5) may be obtained by a relatively efficient elementary 
approach based on the method of weights, but "reduced" with respect to SO (4). Since the 
SO (5) — 7- SO (4) branching rules (25) are known, as is the SO (4) Clebsch-Gordan series, 
the tabulation of weights can be replaced by tabulation of SO (4) irrep labels, which then 
imply all the weights contained within these irreps. To decompose the SO (5) outer product 
(RiSi) ® {R2S2), first the S0(4) irreps in the branchings (RiSi) (XiYi) and {R2S2) 
{X2Y2) are enumerated. Then, for each pair of S0(4) irreps {XiYi) and {X2Y2), the product 
irreps (-^^1^1) (-^^2^2) ^ i^^) ^^^^ enumerated. The aggregate set of these product irreps 
represents the S0(4) content of (RiSi) ® {R2S2). The S0(5) content can now be extracted. 
Namely, the highest weight S0(4) label in the set gives the highest weight S0(5) irrep 
contained in (i^i^i) {R2S2) [which will, incidentally, simply be the sum of (RiSi) and 
{R2S2) as weights]. The SO (4) content of this SO (5) irrep can now be deleted from the 
set, after which the next highest remaining SO (4) label gives the next highest weight SO (5) 
irrep in (i^iS'i) (8) {R2S2), etc. The process is repeated until the set of S0(4) irreps has been 
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(a) 



My 
A 



-0 



(b) 



My 
i 



(c) 



My 



FIG. 3: Branching diagrams for S0(5) D S0(4), shown for (a) the symmetric irrcp (22), (b) the 
antisymmetric irrep (20), and (c) a representative generic irrep (||)- The open circle indicates the 
highest weight for the S0(5) irrep, and the sohd dots indicate S0(4) highest weights, according to 
branching rule (25). The shaded rectangles are the boundaries of the weight sets for these S0(4) 
irreps. The dashed rectangle in panel (c) is the boundary of the S0(4) highest weight set, as 
discussed in the text. 



exhausted. The Clebsch-Gordan series for SO (5) may also be obtained by group character 
methods (see Refs. 18,49). 

The remaining ingredients needed for application of Racah's method are the S0(4)- 

reduced matrix elements of the "additional" generators of SO (5) not contained in SO (4), 
z.e., r(i/2V2). These were obtained in closed form, by solving certain recurrence relations 
obtained from the commutators of the algebra, by Hecht^^ and by Kemmer, Pursey, and 
Williams, as 



{RS) 

(x + iy + i) 



{RS) 
{X + \Y-\) 



T^22> 



T^2 2> 



(RS)\ 
(XY)/ 



{RS)\ 
{XY)/ 



[{R + S - X - Y){R + S + X - 
X {-R + S + X + Y + 1){R 



Y + 3) 

■5 + X + F + 2)]'/' 



2 ' 2 



2X 

[{R + S - X + Y + 1){R + S + X - Y + 2) 
x{R-S-X + Y){R-S + X-Y + l)] 



1/2 



2X + ^Y 



These expressions arc appropriate to the Wigner-Eckart theorem normalization defined in (6) 
and the normalization of T'^^^^^^^^ defined by (24). The remaining matrix elements, connect- 
ing {XY) with (X — iy + i) or (X — — |), follow from these by the self-adjoint property 
of the generators,^'' as 



(RS) 
{XY) 



T^22> 



{RS) 
{X'Y') 



X'Y' ^ ^x-X'+Y-Y' 

\{X'Y') 



{RS) 



T^22> 



{RS) 
{XY) 



(27) 



C. Calculation of SO(5) D SO(4) coupling coefficients 



Consider calculation of the set of S0(5) D S0(4) reduced coupling coefficients for {RiSi)^ 
{R2S2) {RS). It is now straightforward to construct the terms appearing in the Racah 
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condition (11), by use of the chain (I') branching rules, SO (4) Clebsch-Gordan series, S0(4) 
Wigner calculus, and S0(4)-reduced matrix elements of T*^^/^^/^\ compiled in the preceding 
section. For S0(5) D S0(4), the relation (11) may be written 



(RS) 
(XY) 



T^2 2> 



(RS) 

(XT'; 



(XiY,) 



{R2S2] 

iX2Y2] 



(RS) 
(XY) 



= J2 mX,Y,){X,Y,y,{XYm{X[Y;){X,Y,y,{X'Y')] 



\X2Y2 

^22' 



{X[Yi) 
{XY) 



{X'Y') 
(XiY,) 



\X,Y,) {X',Y^) {X'Y') 
{\\) {XY) {X2Y2) 



{RiSi) 
(XiY,) 

R2S2] 
X2Y2] 



T^22> 



T^2 2> 



{RiSi) 
{X'.YD 

{R2S2 
{X'^Yi 



{RiSi) 

(x'.y;) 

RiSi 
XiYi 



{R2S2) 
{X2Y2) 

{R2S2 
{X'.Yi 



{RS) 
{X'Y') 

{RS) 
{X'Y') 



(28) 

The system of equations for the coupling coefficients is assembled following the approach 
of Sec. lie. A different condition is obtained from (28) for each quadruplet of S0(4) 
irreps (XiFi), {X2Y2), {XY), and {X'Y'), chosen from the branchings {R^Si) (XiFi), 
{R2S2) ^ {X2Y2), {RS) {XY), and {RS) {X'Y') [for couphngs involving the identity 
irrep (00), see endnote 34]. A nonnull condition on the coupling coefficients is obtained only 
if {XiYi) ® {X2Y2) {XY) and {X'Y') ® {\\) ^ {XY). Since aU quantities involved in 
the conditions (28) are known exactly, in the form of square roots of rational numbers, and 
since Euler row reduction can be carried out in exact {i.e., symbolic) arithmetic, all coupling 
coefficients can be obtained exactly through the present process, again as (signed) square 
roots of rational numbers. 

Two concrete numerical examples are provided as illustrations of the method in Ap- 
pendix B. A simple low-dimensional example is provided by the coupling (||)(8)(|0) — >■ (|0), 
which involves only a 4 x 4 coefficient matrix, and an example involving an outer multiplicity 
is provided by the coupling (10) ® (l|) — )■ (l|), for which the coefficient matrix has dimen- 
sions 36 X 18. (In the canonical labeling scheme, these examples are [10] ® [||] [||] and 
[11] ® [M] respectively.) 



D. Noncanonical chains 

In the previous section, it was seen how the coupling coefficients for SO (5) reduced with 
respect to the canonical chain may be evaluated by Racah's method, as formulated in Sec. II. 
The necessary ingredients take on a particularly simple form for the canonical chain, in that 
the reduced matrix elements of the generators are given by closed form expressions (26). 
However, the matrix elements of the generators reduced with respect to the noncanonical 
isospin subalgebra [chain (II)] and physical angular momentum algebra [chain (III)] are also 
known. Recurrence relations for the reduced matrix elements have been obtained from vector 
coherent state reahzations,^^'^° or an elementary construction has also been demonstrated for 
chain (11).^^ Hence, Racah's method may be applied directly to the calculation of chain (II) 
and chain (III) reduced coupling coefficients. 

Alternatively, once coupling coefficients for the canonical chain have been obtained, the 
coupling coefficients for the noncanonical chains can readily be deduced by a unitary trans- 
formation. For instance, the couphng coefficients of SU(3) reduced with respect to its 
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physical angular momentum subalgebra are conventionally obtained from the canonical 
SU(3) D U(l) SU(2) coupling coefficients^"*^ '^^ through such a process. The transfor- 
mation brackets between basis states reducing the canonical and noncanonical chains are 
only known in closed form for a few special cases. ^ However, they can be obtained in a 
straightforward fashion, either (1) by diagonalizing the appropriate Casimir operator, i.e., 
or L^, in the canonical basis, as in Ref. 15, or (2) by a combination of laddering and 
orthogonalization operations. The procedure for transformation of coupling coefficients to 
either of the noncanonical chains is discussed in Appendix C. 

IV. CONCLUSION 

It has been shown that Racah's method of infinitesimal generators can be systematically 
generalized to the calculation of reduced coupling coefficients for an arbitrary subalgebra 
chain, provided the matrix elements of the generators (reduced with respect to the lower al- 
gebra) and the recoupling coefficients of the lower algebra are known. For the algebra SO (5), 
the problem of calculating coupling coefficients for generic irreps, reduced with respect to 
the canonical or noncanonical chains, is thereby completely resolved. 

The specific example of SO (5) coupling coefficients may be considered as a prototype for 
the systematic calculation of coupling coefficients for other higher algebras. For instance, 
the computational machinery for SU(3) is well estabhshed^^'^^ and may therefore be used 
as the starting point for calculation of Sp(6) D U(3) reduced coupling coefficients, for the 
fermion dynamical symmetry model, or Sp(6,]R) D U(3) reduced coupling coefficients, for 
the symplectic shell model. The requisite generator matrix elements for Sp(6) and Sp(6, M) 
may be calculated from vector coherent state realizations. ^^'^^ The Sp(6, M) D U(3) coupling 
coefficients are required, for instance, if large-scale calculations are to be carried out in the 
ab initio symplectic scheme of Dytrych et al.^"^'^^ 
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Appendix A: General relations in the presence of outer multiplicities for H 

The derivation of Racah's method in terms of reduced quantities, as given in Sec. II, can 
readily be generalized to the case in which the subalgebra H has outer multiplicities, i.e., its 
Kronecker product is not simply reducible. In this appendix, the necessary generalizations of 
the algebraic relations {e.g., Ref. 18) entering into the derivation of Sec. II are summarized, 
and the fundamental relation (11) for Racah's method is extended to incorporate outer 
multiplicities of H. 
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The general form of Racah's factorization lemma ior G D H is 



oiAi a2A2 
Ai A2 



E 



Ai A2 
Ai A2 



aA 



Tl 12 

aiAi a2A2 



aA 



(AI) 



where a is the multiplicity index for the coupling Ai A2 — > A. The reduced coupling 
coefficients satisfy orthonormality relations 



E 



aiAia2A2 



and 



E 
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aiAi a2A2 
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a'A 



aA 



(A2) 



(^(aiAi)(a;A;)5(a2A2)(a^A^)'5 ITIT' ) 



(A3) 



for any irrep A such that F — > A and, in the second relation, any values of the multiplicity 
indices o and o' . which resolve Ai ® A2 — A and A'^ ® A2 A, respectively. 
For the subalgebra if, the Wigner-Eckart theorem is of the form 




E 



A At 
A Ay 



aA' 
A' 



F' 

a'A' 



F 

aA 



The recoupling coefficients (unitary 6-A symbols) of H are 



Ai A2 0-12A12 

A3 A (723 A23 



E 
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A12A13 



Ai A2 
Ai A2 



(7"l2Ai2 
A12 



A12 A3 
A12 A3 



aA 
A 



A2 A3 
A2 A3 



'^23^-23 
A23 



Ai A23 
Ai A23 



(A4) 



a'A 
A 



(A5) 



This represents the transformation bracket between basis states in the coupling schemes 
[(AiA2)^i2Ai2^g](TA [Ai(A2A3)'"23A23]^'A Under interchange of the ffist and second irreps, 

the coupling coefficients may be expected to satisfy a symmetry relation of the form [e.g., 
for SU(3), see Refs. 51,59] 



A2 Ai 
A2 Ai 



t)-E*-(A^A.^A)(l'li?)- 



(A6) 



Following the same arguments as in Sec. II, the reduced coupling coefficients ior G D H 
are found to satisfy a homogeneous system of equations of the form 
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(A7) 



Tia'SS' 
X 



A2 a; 6' A' 

At a TiAi 



\ aiA] 

TO \ 



a'lA'i 



Ti F2 
a'^A'^ a2A2 



a'A' 



+E 



72 it' 



Ai 




a'A'" 




1 






^ Ti F2 


At 




T2A2 


{ 


ya2A2 




all',) J 


1 aiAi a'2A2 



a'A' 



18 



Coupling cocfBcicnt 

{X-^\\){X.2Y.2){XY) 
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o ^1=^ o 
O '-lies o 



o 
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- 



FIG. 4: Coefficient matrix for the linear, homogenous system of equations determining the S0(5) D 
S0(4) couphng coefficients for {\\) ® (^0) ^ (^0). 



A different equation is obtained for each choice of (aiAia2A2Aa'A') and indices a and r, 
where a resolves Ai (8) A2 — >■ A and r resolves A' A^ — >■ A. 



Appendix B: Numerical examples for the calculation of SO(5) D SO(4) reduced 

coupling coefficients 

For a simple numerical example of the calculation of SO (5) D SO (4) reduced coupling 
coefficients, according to the methods of Sec. Ill C, consider the coupling (||) ® (|0) — )■ (|0). 

The coupling coefficients to be obtained are (^^(qJ)^^ (Jv2) I (Jv2))' ((I/2I/2) (I/2I!) I (iv2))' 

{^^'m^ (1/2 S) I (1/2 S)' and ([1^2 1/2) (iv2) I (1/2?))- The system of four equations in four 
unknowns obtained from (28) has the coefficient matrix shown in Fig. 4. This matrix is of 
rank 3, admitting the null vector 

[-\ -1 I i 1] > 

shown as a row vector for easier comparison with the coupling coefficient labels across the top 
of Fig. 4. It remains to obtain the proper normalization (and phase, if a phase convention is 
to be enforced). The vertical bar demarcates coupling coefficients involving the same product 
irrep {XY), (0|) for the first two coefficients and (|0) for the remaining two coefficients, 
that is, the coefficients which appear together in the normalization sum (15). The squared 
norm, calculated using either pair of coefficients, is N"^ = |. Hence, the normalized coupling 
coefficients are given by 



- — - 

5 V 5 



- - 

5 V 5 



which incidentally also conforms to the generalized Condon-Shortly phase convention. 
One of the simplest cases in which an outer multiplicity occurs is in the coupling (10) ® 

(l|) — 7> (l|), which has multiplicity D = 2. There are 18 coupling coefficients for these 
irreps, related by a system of 36 equations. The first few rows of the coefficient matrix are 
shown in Fig. 5(a). The matrix is of rank 16, admitting two null vectors, given in Fig. 5(b). 
Again, these are shown as row vectors, and groups of coefficients sharing the same (XY) 
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FIG. 5: (a) Coefficient matrix for the linear, homogenous system of equations determining the S0(5) D S0(4) couphng coefficients for 
(10) (8) (1^) The matrix has dimension 36 x 18. Only the first few rows arc shown, (b) Basis vectors for the null space [D = 2). 

(c) Orthonormal coupling coefficient vectors, consisting of the S0(5) D S0(4) reduced coupling coefficients for p = 1 (upper row) and p = 2 
(lower row). 
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are delimited by vertical bars. The inner product matrix (16), which may be obtained using 
any of these four groups of coefficients, is 




The orthonormal null vectors obtained by the Gram-Schmidt procedure with respect to M. 
are then given in Fig. 5(c). The entries of the upper and lower rows may be taken as the 
SO (5) D SO (4) reduced coupling coefficients for p = 1 and 2, respectively. 



Appendix C: Transformation of coupling coefficients 

In this appendix, a general algorithm is outlined for determination of the transformation 
brackets between canonical and UAr(l) <E) S0r(3) [chain (II)] bases of an S0(5) irrep. These 
are the necessary ingredients for deducing the chain (II) reduced coupling coefficients from 
the canonical reduced coupling coefficients. The analogous procedure for transformation to 
the S0l(3) basis [chain (III)] is also briefly considered. 

First, let us review the relevant properties of the proton-neutron quasispin realization of 
SO (5) and the chain (II) basis and branching rules. For protons and neutrons occupying a 
single level of angular momentum j (degeneracy 2j + 1), consider the proton pair quasispin 
operators (X+, Xq, and X_), neutron pair quasispin operators (K|_, Yq, and Y_), proton- 
neutron pair quasispin operators {S^, Sq, and S^), and isospin operators (r+, Tq, and T_). 
Each set spans an angular momentum algebra, which we denote by S0x(3), S0y(3), S05(3), 
or S0t(3), respectively. Explicitly, 



X+ = l{al-al) Xo = -j{al ■ ap + ttp ■ al) X_ = -i(ap ■ a^) 

^+ = 1(4-4) ^0 = -i(4- + 4) Y_ = -l{an-an) 

S+ = ^(Op ■ a.„ + fl^ • Op) 5*0 = ^(O-p ' 5,p -\- dn ' (^n) ^- ~ f (^n ■ Clp + ftp • Cln) 

T+ = -(a],-a„) To = -|(4 • Op - 4 • «n) T_ = -{al-ap), 



where = (-)^-"^a^l, and we define A^^ ■ B^^^ = j{A x B)f = Y.m ^mBm for half-integer 
Thus, for instance, X+ = | Em(-)''""'4,m4 = I Em(4,m«P,m - ap,m4,m)> and 
X_ = I '^^{—y~'^ap,-mO'p,m- There are only ten independent operators, since 5*0 = Xq + Yq 
and Tq = Xo — Yq. The operators defined in (CI) obey the commutation relations of the 
S0(5) generators of Sec. IIIB, with the identifications 

S+ = -2T++ S'_ = 2T__ T+ = 2r+_ r_ = 2T_+, (C2) 

and thus span an SO (5) algebra. 

The generators T+, Tq. T_, and 5*0 span the Uiv(l)®S0T'(3) subalgebra of S0(5), as shown 
in Fig. 2(c). Here Ujv(l) is the one-dimensional algebra of Sq, so denoted in recognition of 
the relation of this operator to the total proton-neutron number operator = Np + Nn. The 
natural weights for chain (II), Ms and Mt, are related to the chain (F) weights by Ms — 
Mx + My and Mt — Mx — My.^^ For the proton-neutron quasispin realization of S0(5) 
defined in (CI), the weights are simply related to proton and neutron occupation numbers 



21 




FIG. 6: Weight diagrams for S0(5) irreps, illustrating the decomposition into UAr(l) S0t(3) 
irreps, i.e., the chain (II) branching. The area of each dot indicates the multiplicity of the weight 
point, and the rectangles indicate grouping of weight points into isospin multiplets. (a) The general 
characteristics, in the presence of branching multiplicites, are illustrated with the irrep (||). The 
branching is obtained by decomposing the weights at a given value of Ms into multiplets, e.g., 
for Ms = +2, in this example, T = 1^, 2^, 3^, 4, and 5, where the exponents indicate branching 
multiplicities. The mathematical labeling schemes for the S0(5) irrep are based on the highest 
weight point (circle), but the {v, t) labeling scheme for pairing applications is based on the quantum 
numbers of the "pair vacuum" UAr(l) ® S0r(3) irrep (dashed box). The arrows indicate the 
action of the isospin laddering operators T± considered in the algorithms presented in the text, 
(b) Decomposition of the low-dimensional symmetric irrep (11). (c) Decomposition for the low- 
dimensional antisymmetric irrep (10), which is the adjoint (or generator) irrep. (d) Decomposition 
of the low-dimensional generic irrep (l|), illustrating isospin multiplets of half-integer isospin. 

by Mx = liNp - Q) and My = i(A^„ - Q), or Ms = |(iV - 2n) and Mr = \{Np - iV„), 
where Q = \{2j + 1) is the half-degeneracy. 

The S0(5) irrep {RS) may be decomposed into Ujv(l) ® S0t(3) irreps labeled by Ms 
and T, as indicated in (17). By inspection of the weight diagram [Fig. 6(a)] and considering 
the relation between canonical weights (horizontal and vertical axes) and chain (II) weights 
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(diagonal axes), it can be seen that the labels are taken from the possible values Mg — 
-{R + S), R + S-1, R + S andT = 01 |, . . ., R + S-1, R + S. For a given value of Ms, 
the highest weight Mt and thus highest isospin is given by Tma^{RS; Ms) = R+S for \Ms\ < 
R-S, and T^^,,{RS; Mg) = 2R-\Ms\ for \Ms\ > R-S. Note that one of the U7v(1)®SOt(3) 
irreps, at the bottom left of the weight diagram, is a "pair vacuum", annihilated by X_, Y_, 
and S'_ [dashed box in Fig. 6(a)]. In applications to proton-neutron pairing, S0(5) irreps 
are conventionally labeled not by the usual mathematical labels (Table II), but rather by 
the seniority v and reduced isospin t^^ These are the occupation number and isospin of the 
pair vacuum, i.e., Msyac = ^{v — 2Q) = —{R + S) and Tyac = t — R — S. 

The isospin content may be seen by decomposing the weights Mt for each given value 
of Ms, i.e., occuring on single diagonal of the weight diagram, into isospin multiplets [solid 
boxes in Fig. 6(a)]. A given pair of U7v(l) ® S0t(3) labels may occur more than once within 
an S0(5) irrep and is therefore labeled by a multiplicity index k — 1, 2, . . ., mult(i?5'; MsT). 
The algorithm for constructing the branching S0(5) — )> U7v(l) ® S0t(3) is derived in, e.g., 
Ref. 8. Simple counting arguments then give a closed form multiplicity formula 

^.mS:MsT) = yf'^'''^^ ^ ^^^-^ (C3) 

^ ^ \f{R + S-T,R- S;Ms) T > R- S, ^ ' 

where f{u, v; w) = Lmin[-u, ^{u + v + w)]\ — [max[0, ^{u — v + w)]] + 1.^^ This relation 
fully defines the branching rule for chain (II). The branching into UAr(l) ® 80^(3) irreps 
is depicted in Fig. 6, for example irreps of S0(5): symmetric [Fig. 6(b)], antisymmetric 
[Fig. 6(c)], and generic [Fig. 6(d)]. 

Basis states for chain (II) involve only a linear combination of canonical chain (I') basis 
states at the same point in weight space. This point may be labeled interchangably either 
by Mx and My or by Ms and Mt. Thus, 



{RS)\ ^ I (RS) {RS)\ (RS) 
I (j^AMxMy Mt I MxM, 



where it is to be understood that Mx = ^i^Ms + Mt) and My = ^{Ms - Mt). The SO (4) 
irreps (XY) of the canonical basis states contributing to a given chain (II) basis state are 
constrained by the branching condition (RS) {XY) and by the usual angular momentum 
projection rules (|Mx| < X and |My| < Y). 

The transformation brackets in (C4) are only known in closed form for a restricted set 
of cases involving low branching multiplicites.^ However, the transformation brackets can 
be systematically calculated for any SO (5) irrep, regardless of multiplicity, through two 
possible procedures, outlined here, involving isospin laddering operations and either matrix 
diagonalization or orthogonalization. 

The first method relies upon the construction of chain (II) basis states as eigenstates of T^. 
Suppose is written as its matrix realization in the (I') basis. The transformation brackets 
appearing in (C4) are the coefficients for decomposition of the eigenstate with respect to 
the chain (I') basis. Therefore, they must constitute the entries of an eigenvector of the 
matrix, of eigenvalue T{T + 1), and may be computed by diagonalization of this matrix. 
Since each chain (II) basis vector involves only a single weight point, the diagonalization can 
be carried out separately on subspaces corresponding to different weight points, in which 
case the eigenvalue T{T + 1) occurs with degeneracy given by mult(it!5'; MsT) from (C3). 
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Note that the matrix reaUzation of at weight point {MxMy) is readily obtained from 
the known matrix elements (26) of the SO (5) generators. From (C2), = (Xq — Yq)"^ + 
2(r+_T_+ + T__|_T+_). This operator can be rccxpressed in terms of spherical bitensor 
coupled products of the generator T^^/^ ^Z^) = (Xq - Fq)^ + 2[(T x T)^^^^ - (T x T)J,J°^]. 

The requisite matrix elements are therefore 



(RS) 
{X'Y') 
.MxMy 



(RS) 
AixI^Yi 



{Mx - My) 



2( (XY) {11] 
MxMy 00 



X 



(RS) 
{X'Y') 



(T X T) 



(11) 



(RS)\ 
{XY)/ 



{X'Y') 
MxMy 



2/ (RS) 

Z{ (XT') 



{T X T) 



(00) 



(|f>). (C5) 



These may be evaluated using Racah's reduction formula/^ as naturally extended to spher- 
ical bitensors, i.e., S0(4) tensors, giving 



{RS) 
{X'Y') 



{T X T) 



(11) 



{RS) 
{XY) 



E 

{X"Y") 



X 



11) 



^22' \22! 

{XY) {X'Y') {X"Y") 



{RS) 
{X'Y') 



T 



{RS) \ / {RS) 
{X"Y")/\{X"Y") 



{RS)\ 
{XY)/ 



(C6) 



and 



{RS) 
{X'Y') 



{T X T) 



(00) 



{RS) 
{XY) 



(XY){X'Y') 



{X"Y") 



/{RS) 
\{XY) 



T 



{RS) 
{X"Y"] 



(C7) 



These expressions involve only SO (4) coupling coefficients (22), SO (4) recoupling coeffi- 
cients (23), and SO (5) D SO (4) reduced matrix elements (26), all of which are readily 
calculated. 

However, it is important to note that the transformation brackets must be consistent 
among weight points within an isospin multiplet, since these are connected by the isospin 

laddering operators T± [Fig. 6(a)]. The chain (II) basis states of a given Ms, k, and T 
but different M^ must be related under laddering by T± with the correct (positive) phase. 
Furthermore, in the presence of branching multiplicities, the choice of basis vectors (reso- 
lution of the multiplicity) must be consistent between weight points, i.e., laddering should 
not connect different k values at adjacent weight points. The relation among transformation 
brackets at adjacent weight points along a diagonal of given Ms, obtained by comparing the 
action of T± on the chain (I') and chain (II) states, is 



[{T ^ Mt){T ± Mt + 1] 



-I 



11/2/ (^y,) 

\m'xM^ 



{RS)\ 
MskT 
M^ 



I (RS) 

= 2 V ( {X'Y') 
^/\M'^M^ 



T. 



{RS) 
{XY) 
MxMy I 



{RS) 
{XY) 
MxMy 



{RS)^ 
MskT 
Mt 



(C8) 



where Mjf = Mx ± |, M(. = My ^ |, and M^ = Mt ± 1, and where it is understood that 
Mx = |(M5 + Mt), My = ^{Ms - Mt), M'^ = l{Ms + M^), and M(. = ^{Ms - M^). 
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At most four SO (4) irreps (XY) contribute to the sum, under the triangularity condition 
{XY) ® (||) — )■ {X'Y'). The matrix element of T-t^ can be calculated from (6) and (26). 

If were diagonalized independently at each weight point, this would not guarantee 
consistency between weight points under laddering. As usual in diagonalization problems: 
(1) Eigenvectors, even in the absence of degenerate eigenvalues, are only defined to within 
an arbitrary phase (sign). (2) In the presence of degenerate eigenvalues (here, multiple 
occurence of the same isospin at given Ms), the choice of orthogonal basis vectors for the 
degenerate eigenspace is arbitrary. Diagonalization of must therefore be augmented by 
further conditions. 

A consistent prescription for the transformation brackets is provided by diagonalizing 
only once, for each Ms, at the most central weight point on the diagonal [Fig. 6(a)], either 
Mt = for integer isospin [Fig. 6(b,c)] or Mt = +| for half- integer isospin [Fig. 6(d)]. The 
remaining transformation brackets are obtained by laddering outward, to more peripheral 
weight points of larger \Mt\ [the "+" arrow in Fig. 6(a)]. If the transformation brackets (C4) 
are considered as entries of numerical eigenvectors, then the laddering operation (C8) con- 
sists of multiplication by the (generally nonsquare and sparse) matrix realization of T-j-^ 
between weight points. For a vector of isospin T, laddering past the weight point Mt = T 
gives a null result. The process terminates at Mt — Tma.x{RS; Ms), or Mt — —Traax{RS; Ms) 
for laddering towards negative Mt- Although the laddering procedure provides a consistent 
set of coupling coefficients, the resolution of multiplicities arising from the diagonalization 
(at Mt = or +i) remains arbitrary and thus still does not provide a unique and repro- 
ducible prescription. Uniqueness (to within phase) can be obtained by further diagonalizing 
a "second" operator, as detailed in Ref. 9.^^ 

The second, alternative approach to constructing the transformation brackets does not in- 
volve explicitly diagonalizing T^. Rather, it is based on laddering inward along a diagonal of 
constant Ms, from the most peripheral weight point [the "— " arrow in Fig. 6(a)] in conjunc- 
tion with orthogonahzation at each weight point. The weight point Mt — Tmax{RS; Ms) 
contains only a single basis state, which therefore is also the T — Tmax basis vector for 
chain (II), to within sign. We are free to choose this sign e.g., as always positive. Lad- 
dering inward to Mt = T^ax — 1 yields the T = Tmax basis vector at this weight point. If 
the number of basis states {i.e., the dimension of the subspace) at this point is larger than 
one, then the remaining orthogonal vector (or vectors) needed to span the space must be 
the chain (II) basis vector (or vectors) of isospin T = Tmax — 1- Ths degeneracy will be 
mult(i?S'; M5, Tmax — !)• The orthogonal basis of good isospin may therefore be found by 
Gram-Schmidt orthogonahzation of a complete but nonorthogonal basis, starting from the 
known T — Tmax basis vector, which is supplemented by mult (i^S*; M5, Tmax — 1) further 
linearly independent vectors as needed to span the subspace. A unique set of transformation 
brackets, both in terms of phases and resolution of multiplicities, is obtained if a well-defined 
prescription is used for specifying these further independent vectors. For instance, one might 
use basis states taken from chain (I'), in order of increasing weight for the S0(4) label (XY), 
starting with the lowest-weight S0(4) label available at the weight point. (Alternatively, 
uniqueness can be enforced by diagonalizing a "second" operator within the T = T^ax — 1 
space, supplemented by a phase convention.) The process of laddering followed by orthonor- 
malization must then be repeated for Mt — T^ax — 2, Tmax — 3, etc., until Mt = or -|-| is 
reached. 

For large-scale computations, the choice between the two methods, (1) diagonalization of 
followed by outward laddering or (2) inward laddering alternating with orthogonahzation. 
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will be dictated by considerations of numerical efficiency and accuracy. These will depend 
upon the numerical linear algebra algorithms being used. The latter method provides the 
simplest route to a unique, reproducible set of phases and resolution of the branching mul- 
tiplicity. 

Once the transformation brackets between bases reducing chains (I') and (11) have been 
obtained, the transformation of reduced coupling coefficients follows immediately. The full 
(unreduced) coupling coefficient may be interpreted as the inner product of a coupled state 
with an uncoupled product of two states, each described by (C4). Then, for the SO (5) D 
Uiv(l) <E) S0r(3) reduced couphng coefficient it follows from the factorization lemma and 
orthonormality that 



{R,S,) {R2S2) 
MsiKiTi M52K2T2 



MskT 



SO(5)DUjv(l)(g)SOT(3) 



E 



(XiYi){X2Y2){XY) 
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(RiSi) 
MxiMyi 
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MxiMyi Mx2My2 



iix^^Y 
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{R2S2) 
{X2Y2) 
Mx2My2i 



(RS) 
MskT 
Mt 



(RS) 
(XY) 
MxMyi 



SO(5)d[Ujv(1)®SOt(3)oSO(4)] 



X 





{R2S2) 




(X2Y2) 



SO(5)dSO(4) 



for any value of Mt allowed given isospin T, where again it is to be understood that 
Mx = l{Ms + Mt), My = ^{Ms - Mt), and similarly for Mxi, Myi, Mx2, and My2- 
The sum over S0(4) irrep labels is subject to the usual branching and triangularity con- 
straints on the canonical reduced coupling coefficients. For the resulting chain (11) reduced 
coupling coefficient to be nonvanishing, it must obey the S0(5) — )■ Ujv(l)® 80^(3) branching 
rules (C3), the UAr(l) additivity condition Msi + M52 = Ms, and the 80^(3) triangularity 
condition Ti T2 T. 

In Appendix B, the canonical reduced coupling coefficients for the S0(5) coupling (10) ® 
— > (l|) were calculated [Fig. 5(c)], as a relatively simple numerical example involving 
an outer multiplicity [D = 2). For a concrete illustration of the transformation procedure 
just described, let us consider the transformation of these coefficients into chain (II) reduced 
coupling coefficients. The decomposition of the irrcps (10) and (li) into Uiv(l) ® S0t(3) 
irreps is shown in Fig. 6(c,d). For the transformation brackets, wc follow the second method 
above. For instance, consider the Ms = +| diagonal of the weight diagram for the (l|) 
irrep of S0(5) [Fig. 6(d)]. The T = | seed state at Mt = +| is | = +1 (^D+l-i)' 

where we abbreviate chain (I') basis states as | {XY)MxMy) and chain (II) basis states as 
I MsTMt). Laddering inward to Mr = -F^ yields T = f state 



= +vil (iO)+io) + J\\ (ii)+io). 



Orthogonalization, using | (|0)-|-|0) as the independent Gram-Schmidt basis vector, gives 
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TABLE III: Chain (II) reduced coupling coefficients for the S0(5) couphng (10) (g) (1^) (1^), 
obtained from the canonical couphng coefficients of Fig. 5(c) according to the transformation 
conventions described in the text. For each UAr(l)(8)SOT(3) coupling {Ms iTi)®{Ms2T2) {MgT), 
the coefficients are given for S0(5) outer multiplicity index values p = 1 and 2, respectively. 

Msi Ms 2 Ms Ti T2 T Coefficients M51 Ms 2 Ms Ti T2 T Coefficients 



1 - 



1 
1 
1 
1 
1 
1 
1 
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1 
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J7 

'16 
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T = I State 

I = - Til 

Continued laddering to negative Mt gives the remaining transformation brackets along the 
diagonal. Straightforward apphcation of (C9), using the full set of transformation brackets 
derived in this fashion, then yields the chain (II) reduced coupling coefficients in Table III. 
Note that the coefficients are either symmetric or antisymmetric under the Ujv(l) particle- 
hole conjugation operation Ms —Ms (see Ref. 9). 

The S0l(3) subalgebra of chain (III) is the maximal S0(3) subalgebra, i.e., one which 
is not contained within any larger proper subalgebra of SO (5). This subalgebra is obtained 
by letting^^-^^ 

= -V2X+ - y6T„+ = Xo + 3^0 l'^-I = V2X_ + V6T+_ 

= -^/5Y+ O^^l = V1OT++ = -V3X+ + 2r_+ OS^^ = 3X0 - Yo (CIO) 

Oi^J = V3X_ - 2T+_ O^^l = -\/lOT__ O'lI = V5y_. 
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Then 80^(3) has generators L^Ml-i where, as usual, L^i 

O^^ constitute an octupole tensor with respect to S0£,(3). The commutation relations of 



^^L±. The remaining generators 



the generators are given (most compactly in spherical tensor coupled form^^) by 

[L, L] ^-V2L [L, O] = -2^/30 
[O, Of^ = -2V7L [O, Of^ = ^O. 



(Cll) 



The S0l(3) weight Ml is related to the canonical weights by Ml = Mx + 3My, according 
to (CIO), and thus defines an oblique axis in the weight space. The S0(5) generators are 
shown classified according to this weight in Fig. 2(d). The dashed fines connect the canonical 



generators of Sec. Ill B, which have good Ml but not L, to the finear combinations L^^^ 

and 0\/^, which have both good Ml and good L. 

The labels for basis states reducing the 80^(3) subalgebra are indicated in (17). Basis 
states for chain (III) involve a linear combination of chain (I') basis states of the same Ml- 
The transformation between bases therefore involves a sum not only over {XY), as in (C4), 
but also over distinct weight points (MxMy), 
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(C12) 



subject to the constraint Ml = Mx + 3My. The transformation brackets may be systemat- 
ically evaluated using an adaptation of either of the methods proposed above for Chain (II): 
(1) diagonalization of for the Ml = or +^ subspace of the S0(5) irrep, followed 
by laddering outward to larger-M^ spaces, or (2) laddering inward from the maximal-M^, 
space, in conjunction with orthogonalization within each successive lower- M^ space. For a 
unique resolution of the S0(5) D 80^(3) branching multiplicity, it has been suggested that 
chain (III) basis states be chosen in which the octupole tensor is diagonal, equivalent to 
diagonalizing the Hermitian "second" operator {L x L x O x L)(°^ + {L x O x L x L)(°).^® 
Once transformation brackets have been obtained, the transformation of reduced coupling 
coefficients is given by 
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SO(5)dSO(4) 
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for any value of allowed given angular momentum L. where the summations over Mxi, 
Myi, Mx2, My2, Mx, and My are subject to the constraints Mn — Mxi + 3Myi, Ml2 — 
Mx2 + 3My2, and Ml = Mx + SMy. 
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Eulcr row reduction permits solution for the null vectors in exact symbolic arithmetic, as il- 
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